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r-j I Abstract We discuss on the construction of the Lagrangian for a classical sys- 

P-| tern with mass varying explicitly with position. The conventional Lagrange 

^ equations cannot be directly applied to this system because there is a non- 

conservative force quadratic in the velocity which is associated to the variable 
mass. We construct a Hamiltonian for this system and find the modifications 
which are required on the Euler-Lagrange and Hamilton's equations in order 
to reproduce the appropriate Newton's dynamical law. A canonical transfor- 
mation is found to map the variable mass equations to the ones of a constant 



mass. The time-dependent first integrals of motion are found in terms of the 
factorization of the Hamiltonian. As particular cases, we recover some recent 
results for which the dependence of the mass on the position was already un- 
noticed, and find new solvable potentials of the Poschl- Teller form which seem 
to be new. The latter are associated to either the su(l, 1) or the su{2) Lie 
algebras depending on the sign of the Hamiltonian. 

^ ■ 1 Introduction 

In recent papers [IH1] the factorization of the classical Hamiltonian in terms of two 
functions that, together with the Hamiltonian itself, close a Poisson algebra was dis- 
cussed. This program included mechanical systems with mass explicitly dependent on 
position [U[2] which was extended to the quantum case [UI21E1E]. Yet, the discussion is 
far from being exhausted. The mass-function m(x) gives rise to a 'force quadratic in the 
velocity' which, in turn, leads to a non-linear differential equation of motion. The most 
celebrated example of this kind of equations is due to Mathews and Lakshmanan [7] (see 
also [E1IH])- They analyzed a non-linear oscillator obeying the equation 

(1 + \x 2 )x - Xxx 2 + a 2 x = 0, (1) 

and showed that this last equation is derivable from the Lagrangian 
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In this paper we discuss on the construction of the Lagrangian for position-dependent 
mass systems (Eq. [2] is intended to be a particular case). In Section [2, we depart from 
the Newton's equation of motion for these systems and show that it is possible to get a 
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Lagrangian in the standard form L = 2 m(x) ~ "^(• E )' a ^ the cos ^ °f modifying the Euler- 
Lagrange equations of motion. The modification must include the non-inertial force which 
is quadratic in the velocity. Then, the Hamiltonian is constructed and the canonical 
equations of motion are shown to be consistent with both, the Newtonian and the La- 
grangian formulations already discussed. Although the Lagrangian and Hamiltonian are 
not unique, our standard expressions are found to be in correspondence with the general 
expressions derived by other authors [10J. Some other results already reported in the 
context of constant mass systems in presence of forces quadratic in the velocity turn out 
to be particular cases of our position-dependent mass approach. We extend the results 
on the factorization method [IHU to the case of the variable mass systems obeying the 
new canonical equations of motion. We find that only potentials of the Poschl- Teller 
form allow the construction of time-dependent first integrals if the underlying Poisson 
structures have to close a deformed Poisson algebra. Then, explicit expressions for the 
phase trajectories (x(t),p(t)) described by these variable masses are found. The section 
is closed with the identification of a point transformation leading from the variable mass 
problem to the one of a constant mass. This transformation is shown to be canonical in 
the sense that it leaves invariant the Hamilton's equations of motion here derived. The 
Section [3] is devoted to the analysis of the diverse Poschl- Teller potentials we have at hand 
by supplying specific mass-functions m(x) in our approach. In particular, we use a doubly 
singular mass m ss (x) arising from the analysis of the inversion potential for NH 3 in terms 
of the density operator [11]. The potentials arising from a singular as well as a regular 
masses which we have used in our previous works are also analyzed. Finally, an exponen- 
tial mass-function m e (x), introduced here to recover the Lagrangian forms reported by 
other authors for constant mass systems, gives rise to solvable potentials which, as far as 
we know, have been not reported previously in the literature. The paper is closed with 
some conclusions. 



2 Deformed algebras of a particle with position-dependent 
mass 

2.1 Newtonian framework 

For the sake of motivation consider the one- dimensional motion of a system of mass 
m(x) > 0, explicitly dependent on the position x, acted by a force F which in general 
depends on position x, velocity x, and time t (the case of n generalized coordinates is 
straightforward). The Newton's equation of motion is 

F(x, x; t) — -j- = m'(x)x 2 + m(x)x, (3) 
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where p = m(x)x is the linear momentum. Hereafter g and g' will stand for time and 
position derivatives of g respectively. To get fl3]) we have assumed a null velocity v m for 
the accreted or emitted mass relative to the system (otherwise, one of the velocities in 
the quadratic term m\x)x 2 must be substituted by the displacement v m + x). Equation 
02]) has the general form 

F(x, x; t) — a(x; t)x + /3(x, x; t), (4) 
and can be decoupled in the nonlinear system 

x = v, v — — ^-r- \F(x, v; t) — m'(x)v 2 ] . (5) 
m[x) 

Equation (jSJ) defines a path on the plane (x, v) in terms of the parameter t (see, e.g. [32] )• 
The latter expressions acquire the standard form whenever ml = 0: 

x — v, v = — F (x,v;t). (6) 
m 

Here the constant m > corresponds to the mass of a system subjected to the action of 
the force F . It is possible to make the equivalence between the two systems, (jSJ) and (jSJ), 
in such a way that the variable mass m(x) describes the same path as a given constant 
mass m . For instance, consider the mass 

m(x) = — m (l — kVo(x)) (7) 

with k a constant, and Vq(x) a differentiate function. Let E = 1/k be the total energy 
of m , then Eqs. §5§ and §6§ lead to equivalent paths provided that 

2E p 2 dV»(x) (g) 



m v 2 dx 

Notice in particular that the variable mass (J7J) describes, under its own inertia (i.e., in 
absence of the external field F), exactly the same path as the constant mass mo under 
the action of the potential function 2Vq(x). This example shows that the mechanical 
energy Eq of the properly chosen constant mass system can be used to parameterize the 
path described by the variable mass m(x) on the (x, v )-plane (regardless the mechanical 
energy is not conserved in the latter system). However, we have to emphasize that this 
equivalence holds for the paths only, not for the motion with respect to time |13j . 

To get some insight on the forces involved in the Newton's equation let us rewrite 
this in the standard form 

m(x)x = F net (x, x; t) = F(x, x; t) — m\x)x 2 . (9) 

Since x 2 > 0, the system is accelerated if the rate m'(x) is negative, and decelerated if 
this is positive. Thus, the term quadratic in the velocity —m'(x)x 2 corresponds to the 
thrust of the system, and indicates how this alters the velocity x. In this way, the net 
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force F net on a particle suffering a spatial variation of its mass results as a combination 
of the external force F and the thrust — m'(x)x 2 . 

Mechanical systems of variable mass are mostly studied in the cases where the mass 
m is an explicit function of time m(t). Immediate applications involve rockets and jet 
engines analyzed in terms of Newtonian theory. The main difficulty is that not all the 
involved forces are derivable from an ordinary potential function or even from a gener- 
alized potential. In many situations of interest, the action of such forces produces the 
non-conservation of the mechanical energy (see for instance the analysis of the raindrop 
problem in [14J). Nonetheless, the dynamics of time- dependent mass systems can be 
studied in the Lagrangian approach. Indeed, for these systems in one degree of freedom 
and forces independent on velocity, Darboux showed that it is always possible to con- 
struct a Lagrangian [15] (see also [IS]). The equation of motion for some of these systems 
is connected to that of a pendulum whose length £ varies with time in the small-angle 
approximation £9 + 2£6 + g9 = (0 < 9 << 1). For instance, a time-dependent mass 
(undamped) oscillator obeys the equation mx + mx + kx = 0. Here, the counterpart of 9 
is the displacement x for the oscillator and the thrust —fax simulates a damping linear in 
the velocity. It is well known that the solutions of this last equation are either decreasing 
or increasing oscillatory functions, depending on the sign of m. Interestingly, a simple 
experimental setup can be used to show that the standard expression for the oscillator's 
frequency u = y/k/m remains valid if the mass is time-dependent p2]. A useful result 
since the mathematical tools used to analyze the constant mass oscillators can also be 
applied if the mass depends on time. Our approach considers a mass which is an explicit 
function of position m(x), so that implicit dependence on time is assumed m(x(t)). Ac- 
cording to ([3]), the equation of motion of a position-dependent mass oscillator is nonlinear 
m(x)x + m'(x)x 2 + kx = (compare with Eq. [1]). In contraposition to the time-dependent 
mass case, the thrust —m'(x)x 2 is here quadratic in the velocity. In general, it is not 
evident how a set of canonical variables can be found such that the Newton's equation 
(J3J), even in the latter oscillator's case, is fulfilled. 

In the sequel, we shall establish a mechanism to determine the phase trajectories 
(x(t),p(t)) associated to the Newton's equation of motion (j3j) for interactions F derivable 
from a properly chosen potential function. We will focus on the invariants which arise after 
factorizing the related Hamiltonian. With this aim, we first analyze the problem in the 
Lagrangian and Hamiltonian frameworks to determine whether the standard expression 
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of the Hamiltonian H = + V remains valid if the mass depends on position. 
2.2 Canonical framework 

Departing from the Newton's equation of motion ([3]) for a force independent on velocity, 
and applying the D'Alembert's Principle we arrive at 




(10) 
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with T and R the kinetic energy and reacting thrust, respectively given by 

T:=-m(x)x 2 , R(x, x; t) := — -m'(x)x 2 . (11) 

When the external force is derivable from a scalar potential function F = — and this 
last does not depend on either velocity or time, equation ({TO]) becomes 

d fdL\ dL ~ 

This last is the Lagrangian form of the Newton's equation of motion (J3]). To verify that 
the Lagrange equation (fT2l) encodes the dynamics of the variable mass m(x) we are dealing 
with, consider that the acceleration a = x arising from L satisfies 

Ma = Q, (13) 

with M = the Hessian (matrix) of L with respect to the velocity, and the "Lagrangian 
force" Q defined as 

Q = ^L^ ± ^L + R. (14) 



dx dxdx 

dx 



The force Q includes the applied (physical) force and the fictitious (non-inertial) 

force lr — Air^r = —R- This last shows that the reacting thrust R is a non-inertial force. 

ox axax 

The substitution of L = ^m(x)x 2 — V(x) in f lTSlTHl) reproduces the dynamical law (jS])- 
Moreover M = m(x) > 0, so that the Lagrangian L is regular and equation ffT3"|) can be 
solved for the acceleration a = Q/m(x). The solutions x(t) of the second order differential 
equation 'x = a(x, x; t) can be constructed in terms of the non-linear system ([5j). 

Once we have constructed the Lagrangian L for the posit ion- dependent mass m(x), 
we use the momentum p = ^ = m{x)x to obtain the Hamiltonian H from the Legendre 
transformation 

p 2 

H(x,p; t) = px — L(x, x; t) = - — -r— + V(x). (15) 

2m(x) 

The Hamiltonian's time rate of change 

at m[x) 

with the reacting thrust R expressed in canonical coordinates 

R = R(x,p;t) =R(x,x(x,p;t);t) = -^M (^-^j , (17) 

shows that the value S of % is not independent of t. On the other hand, a simple 
calculation leads to the appropriate canonical form of equation ( IT2l) . One gets 

x = — , p =-— + R. 18 

op ox 
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The differential equations (fT8l) correspond to the law of motion of the mechanical position- 
dependent mass system in phase space. At time t, their solutions are the canonically 
conjugate variables of position and momentum {x,p} = 1, with {, } the Poisson bracket 

r, i df_dg_ _ dj_dg_ 

' dx dp dp dx 

The Hamilton's equations f JT8|) in the Poisson bracket formulation read 

x = {x J U}, p={p,U} + R. (20) 

In general, the explicit form of m(x) and V(x) determines the Hamiltonian's domain of 
definition U(H). A point transformation is then expected to reduce equations (jTSJ) to the 
ones associated with an equivalent constant mass system in the appropriate domain (see 
Section E3D- 

To summarize this section, let us rewrite the Lagrangian and the Hamiltonian associ- 
ated to the Newton's equation of motion ([3]), we have 

1 » 2 

L = -m{x)x 2 -V{x), %= -^— + V{x). (21) 

2 2m(x) 

These results show that the standard expressions for L and % remain valid if the mass is 
an explicit function of the position. Thus, a simple description of the position-dependent 
mass particle interacting with an external environment consists of replacing the (constant) 
mass m by the appropriate function of the position m(x) in L and H. Although the 
Lagrangian and Hamiltonian are not unique, our expressions (1211) are ensured by the 
conditions for the existence of standard Lagrangians for equations with space-dependent 
coefficients discussed in [10]. Let us multiply Eq. (J3J) by l/m(x), after the identification 

m'{x) dlnmjx) 1 dV{x) 

Kx) = — r~r = , c(x) = —-, g{x) = — — , 22 

m[x) dx m(x) dx 

we arrive at 

x + b(x)x 2 + c(x)g(x) = 0. (23) 

So that this last equation admits a Lagrangian description (see Proposition 3 of pi]]), and 
its standard Lagrangian is 

L(x,x) = ^x 2 e Ib(x) - J c(x)g(x)e Ib{S;) dx, (24) 

where the quantity 

h(x) = / b(x)dx = lnm(x) (25) 



reduces fl24l) to our Lagrangian (1211) . This standard form of writing L and H is appropriate 
to recover (as particular cases) some of the Lagrangians and Hamiltonians already reported 
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in a different context by other authors. For example, consider a mass-function m e (x) and 
a potential function V(x) such that 

m e (x) = m e^\ ^ = e «/»^ ) (26) 

ax ax 

with mo and l//c constants expressed in mass and position units respectively. The New- 
ton's equation (J9j, Lagrangian and Hamiltonian (121]) become 



aV(x) 1 2 

m x = -m nx , (27) 

ax 2 



L = -m xV*/ 2 - f dse KS ' 2 "H = + f dse KS ' 2 < ^-. (28) 

2 Jo « s 2m o Jo ds 

These last expressions are in agreement with those reported in [18] for a system of con- 
stant mass m , subjected to a 'force quadratic in the velocity' — |m /tx 2 . An immediate 
generalization, considering now a mass-function m(x) = m exp(|/t(x)), can be put in 
connection to the system of constant mass discussed in [19] . It is remarkable that both of 
such results arise from the dynamics of a particle with mass varying exponentially with 
the position, a situation which seems to be unnoticed in Refs. [181 119]. Moreover, no 
solutions to the corresponding canonical problem ( 120]) are given in neither [18] nor [T9] . 
Next, we shall construct solutions to the system ( [20]) for m(x) an arbitrary, differentiable 
and integrable function of x, and the properly chosen potential V(x). 



2.3 Factorization and deformed algebras 

The dynamical problem (120]) can be studied in two general ways (compare with the quan- 
tum problems studied in [5]). First, given a specific potential V(x) acting on the mass 
m(x), the related phase trajectories are found. Second, given an algebra which rules 
the dynamical law of the mass, the potential and phase trajectories are constructed in a 
purely algebraic form. Our aim here is to follow the second approach. For this we shall 
extend the factorization method discussed in [3JH] to the case of classical systems having 
a position-dependent mass and obeying the dynamical law ([3]). The factorization of the 
Hamiltonian ( 1T5]) leads in a natural form to the identification of a pair of time- dependent 
integrals of motion Q which, in turn, allows to construct the phase trajectories (x(t),p(t)) 
associated to the canonical equations (TT8j) . 

We look for a couple of complex functions A + (x,p;t), A~(x,p;t), and a constant e 
such that the Hamiltonian ( [15]) becomes factorized 

H = A + A- + e = A'A + + e. (29) 

The inspection of these last equations suggest to define A ± as follows 

A ± = Tif(x) J-—^ + g(x)^fU, (30) 
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with f(x) and g(x) functions to be determined. Here, we are considering the possibility 
of bound states (confined motion) for which the Hamiltonian H is of negative sign % < 0. 
Thereby, according to the sign of H, 7 will be either +1 or —1 such that \pYH. is real in 
( ]3"U]h The introduction of f[5U|) into ( 125]) leads to the relationships 

f 2 (x)+ 19 2 (x) = l, (31) 

e = V(x)/ 2 (x). (32) 

Given m(x) and V(x), the functions A^ and the Hamiltonian H induce a Poisson struc- 
ture, 

i{A-,A+} = —L=W(f,g)y/lH, 

\/2m{x) (33) 

{n,A+} = A+{A-,A+}, {n,A-} = -{A-,A + }A-, 

where W(f, g) = fg' — f'g is the Wronskian of / and g, and we have used (132]) . Now, we 
ask the system ( 133]) to close a deformed Poisson algebra by demanding that {A~ , A + } be 
expressed in terms of the powers of \/yH. The easiest way to satisfy this condition is by 
looking for the solutions of equation 



W(f,g) = a-\/2m(x), a = const. (34) 
Using f j3T]) . this last equation is reduced to quadratures 



I , d9 = v/2o%^ [ X J(t)dt, J(x) = (35) 
J a/1 - ig 2 Jc V m ° 

The constants mo and c are expressed in mass and length units respectively. Considering 
the sign of the Hamiltonian H, equation f ]35]) gives 

r sin [V2a 2 m f*J(t)dt] , 7 = 1 (H > 0) 

9(x) = I " (36) 

[ sinh [V2a 2 m f*J{t)dt] , 7 = -1 {% < 0) 

The expression for /(x) is obtained after the substitution of (]36j) into (13B . At this stage 
we realize that the relationships ( 13111321) define the potential V in terms of the ^-function 



V(x) 



{n > 0) 



cos 2 [V2a 2 m J x J(t)dt] 

(37) 

^ cosh 2 [V2^r f*j(t)dt\ {n > 0) 



Thus, the potential allowing equations ( 151]) and (133]) is not arbitrary. As we can see, given 
the mass m(x), a potential of the Poschl- Teller form (13T]) is such that the system (1331 



S 



becomes the deformed Poisson algebra defined by 



i{A~,A + ) = 2«V7?7, (38) 

i {%, A*} = ±2a^iKA ± . (39) 

The latter results are consistent with our factorization approach ( l29tl3T)|) . Indeed, using 
relation ( 139]) . it is easy to verify the following Poisson brackets 

{H, A + A- } = {H, A'A + } = 0. (40) 

Hence, the factorization f l29|) makes sense because the products A + A~ and A~ A + are 
functions of the Hamiltonian %. In particular, up to an arbitrary additive constant, each 
of them can be chosen to be proportional to H. Moreover, using ( 139]) and the Jacobi 
identity 

{H, {A-,A + }} + {A-, {A + , H}} + {A + , {H, A'}} = 0, (41) 
one arrives at the Poisson bracket 

{H, {A~, A + }} = -i2a{A + A-, y^fH} = 0. (42) 

So that {*4r,.4 + } is also a function of H and (138]) makes sense. 

The relevance of equation ( 139]) is clear, it implies both the factorization ( 129]) and the 
Poisson bracket (138]) . This also allows the construction of first integrals of motion for the 
problem we are dealing with. Let us introduce the functions Q = 0±(t)A ± , and calculate 
their total time derivative 



30 

±i2ay/iH6±{t) + 

at 



A ± . (43) 



The ^-functions are determined by canceling the expression in square brackets. For S 
the value of H, taken as a parameter in the context of Section I2.1[ one gets 0±(t) = 
#oe =R2av/ ^' (hereafter 0q = 1). These time-dependent functions also produce a null total 
time derivative (|43p . so that are the following time-dependent integrals of motion 



Q ± = exp [Ti2a y/^S t] A ± . (44) 

Remark that Q + is the complex conjugation of Q~, therefore Q + Q~ = \Q~\ 2 — % — e. 
Let q* 1 be the value of the integral of motion Q , in polar form this can be written as 



q ± = VS-e e ±i *°, O G R. (45) 

Let us stress that £ — e > implies £ > e > for % > 0, and \£\ < \e\ for H < 0. The 
values of A ± can be now obtained from (13"D|) . OS]) and (H5|) : 

T if{x)-=E== + g{x)^£ = V£~^~e exp[± ? (0 o + 2a V^*)]- ( 46 ) 
y 2m{x) 
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The latter equations lead to the phase trajectories (x(t),p(t)) of bound states (confined 
motion) with energy S greater than or equal to the global minimum of the potential. The 
constant </>o must be fixed by the initial conditions. We have 



g(x(t)) = \\—jr cos[0 o + 2a v / 7^t], (47) 



where 



P(*) = ~\l l(t) sm[<f>o + 2ay^et}, (48) 



The explicit form of x(t) is obtained from (1471) and ( 1361) by using the composition [g 1 
9Mt)). 



2.4 Point transformations 

In Section 12.31 we have shown that potentials of the Poschl- Teller form ( 137|) give rise to 
the deformed Poisson algebras ( 13811391) ; these last associated to the dynamics of a particle 
of mass varying with the position. Now, we want to make a point transformation leading 
from the canonical coordinates, x and p, of the variable mass m(x), to the coordinates, 
Q and P, of an equivalent system of constant mass m . As usual, we shall write 

Q = Q(x,p;t) and P = P(x,p;t) (50) 

for the point transformation we look for. Using (fl8l) we arrive at 

+ + * P ={PlVi + B »L + %. (51) 

We want these equations to have the same form as those given in (fl8|) for a function 
K(Q, P; t) to be determined. Thus, we ask for 

• OK . dK n 

Q= 8P' P= -dQ +R '" (52) 

where Rk is the reacting thrust (TPTj) expressed in the new coordinates. For the sake of 
simplicity, let us assume that Q is not an explicit function of neither p nor t, and that P 
is time- independent. We have 

dK dK dP 
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Let K be the Hamiltonian function associated to a system of constant mass mo under the 
action of a given potential V(Q), 



K = J^- + V{Q). (54) 



2m r 



If V(Q) is of the Poschl- Teller form 



cos 2 (a Q) 
e 



V(Q) = 

k cosh 2 (a Q) 

one easily identifies that the point transformations 



(K > 0) 

(55) 

(if < 0) 



Q 



[ J(s)ds and P = -j— (56) 

J J [Xj 



allow the mapping from the variable mass Hamiltonian H to the new one K. It is a matter 
of substitution to verify that these new variables satisfy the equations ( 1521) . explicitly 

. p . dv(g) 

<5 = — , P = + Rk- (57) 

mo aQ 

That is, equation ( 156]) represents the canonical transformation from the dynamics of the 
position-dependent mass Poschl- Teller potential A3 Tj) to the one of the constant mass 
system (j55p . This result is consistent with the phase trajectories 047H48p . for which the 
value £ of H is taken as a parameter. 

3 Position-dependent mass Poschl- Teller potentials 



In the previous sections we have shown that potentials of the Poschl- Teller form fl37|) lead 
to the solving of the canonical equations of motion f|T8|) in terms of the invariants defined 
by the deformed Poisson algebra (I38H39I) . Depending on the sign of the Hamiltonian ~H, 
equation fl3Tl) includes two general forms of these potentials and, as a consequence, we also 
have two different realizations of the algebra (I38ti39p . Namely, for a positive Hamiltonian 
% > 0, the identification 

a* = -A ± , a = -Vn, (58) 
a a 

makes clear that the Poisson structure (I38H39I) associated to the dynamical algebra of the 
trigonometric position-dependent mass Poschl- Teller potential 

v(x) = cos 2 [V2^r ffjttjdt\ (59) 



11 



corresponds to the su(l, 1) algebra 



i {a",a + } = 2a°, z {a , a ± } = ±a ± . (60) 
In a similar form, for a negative Hamiltonian % < 0, the identification 

a ± = -A ± ) a = --V^n, (61) 

a a 

shows that the dynamics of a position-dependent mass particle subjected to the hyperbolic 
Poschl- Teller potential 

v(x) = cosh 2 [V2^r j*j(t)dt] (62) 

is nothing but the su(2) Poisson algebra 

i{ a -,a + } = -2a°, i{a°,a ± } = ±a ± . (63) 

Next, we are going to discuss some specific realizations of these results by supplying 
different forms of the mass-function m(x). In each case, the phase trajectories (x(t),p(t)) 
are explicitly derived. 

3.1 Doubly singular mass-functions 

The mass function 

m ss (x)=m (^ 1 1 ~^Q , A>0, < (64) 

has been used in the study of the inversion potential for NH3 in terms of the density 
theory [11]. This function is singular at the points x = ±1/A, and has no zeros in the 
interval (— 1/A, 1/A). In notation of Ref. [IT] , m = fio should correspond to the reduced 
mass 

3mM 
1X0 = 3m + M' 

with m the mass of each of the three hydrogen atoms, and M the mass of the nitrogen 
atom in the geometry of a rectangular pyramid where M is in the cusp. The mass-function 
(1641 then corresponds to the situation in which the distance between each pair of masses 
m is allowed to change. The parameter 1/A = ro stands for the separation between m 
and M at the planar equilibrium geometry, x is the inversion coordinate and the negative 
parameter /3 is given by /3 = —3m/M. Here, m > and 1/A > are arbitrary constants 
expressed in mass and distance units respectively. The negative, dimensionless parameter 
(3 is also arbitrary. 

To get explicit forms of the solutions we use (l35i) and arrive at 



C J ss {t)dt= f\h^dt= \ r \l l —^rdt = \e (arcsin(Ax),/3 1 /2) 
Jo v m o A Jo v 1 — t A 
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(hereafter we take c = 0), with E(<f),k) the Incomplete Elliptic Integral of the Second 
Kind [20]: 



Potentials (J37J) are in this case given by 



cos 2 [y / 2moyi?(arcsin Ax, /3 1 / 2 )] 



(H > 0) 
(H < 0) 



(65) 



cosh 2 [y^mof E (arcsin Ax, (3 1 / 2 ) 
The domains of definition for these potentials are 

V ss> = {x | -| < V^jE (arcsin Ax, /3 1/2 ) < | } , (66) 

and 'P ss < = (— 1/A, 1/A) respectively. Finally, the corresponding phase trajectories are 
obtained from the expressions 



sin y/2m —E (arcsin Ax(t), P 1 ^ 2 ) = y — — — cos ^2,VSat + </> j , 



P(t ) = ^(S- t)( i-B<MW) sin , 24lat + 

1 - (Ax(t)) 2 cos [y/2m^fE (arcsin Xx(t), /3 1 / 2 )] 



(67) 



and 



sinh y/2m —E (arcsin \x(t), Z? 1 ^ 2 ) = y — — cos ^2V— + 0o 

P(i ) = V2,tf- t )(l-W4)ff sin / 2v ^ tti + ^ 

7 1 - (Ax(t)) 2 cosh [v 7 ^! ^ (arcsin Ax(t), /3V2V 



(68) 



The potentials V ss (x) and phase trajectories (x(t),p(t)) have been depicted in Figures [T] 
and |2] for different values of the parameters. From Figure [1] (left) we realize that only 
confined motion is allowed for the domain T> ss> . The mass m ss (x) is singular at the edges 
of T> ss> , where the particle reaches the turning points and its momentum p is zero. The 
minimal amount of mass is obtained at origin, so that the mass increases in value as 
the particle approaches the turning points and conversely, the particle losses mass as it 
approaches the origin. The momentum p evolves in time in a quite similar manner: For 
bounded energies greater than the global minimum of the potential, the particle acquires 
a momentum p(t) which is greater as the particle approaches one of the turning points. 
Once there, the momentum of the particle changes in sign so that the motion is reverted 
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by a strong acceleration. This relationship between the maxima of the mass and the 
maxima of the momentum is also found in the confined motion for % < (Figure |2l 
left). There, it is also true that the momentum p is zero at the points where the mass is 
divergent. The same can be said for the other masses analyzed in the sequel. Of particular 
interest, the regular mass case (discussed in Section I3.3P is such that the momentum p 
and the mass-function mu(x) are maxima at the origin, and both of them take their lower 
admissible values at the corresponding turning points. 



V(x) V(x) V{x) 




Figure 1: Upper line. Position dependent- mass trigonometric Poschl- Teller potentials 
for the mass-functions m ss (left), m s (center), (right), and the parameters a = 1/2, 
= —1, m = A = 1. Lower line. Some of the related (confined motion) classical 
trajectories in phase space for </> = 0, e = 2 and £ = 2.5 (red curve), £ — 4 (dashed 
curve), £ = 5.5 (blue curve). 



3.2 Singular mass-functions 

Given a mass function m(x), the quantum problem for the potential V(x) can be solved 
by the mapping of the Schrodinger equation of m{x) to the Schrodinger equation of a 
constant mass mo- In the trivial case of a constant mass mo, the new potential V e ff 
(expressed in the new coordinates) is the same as the former one V. In general, the mass 
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function 



m s (x) 



m 



(69) 



(1 + Ax) 2 

has been shown to be the simplest nontrivial case in which V e ff = V [5]. This mass is 
also connected with the revival of wave packets in a position-dependent mass infinite well 
], and was used in our previous analysis of diverse position-dependent mass oscillators 
The function (169]) is singular at the point x = — 1/A and has no zeros in 
[— 1/A, +00). The potentials f )37|) read as 



VJx) 



cos^ [VZfm f ln(l + Ax)] 



cosh 2 [V2m^f hi(l + Ax)] 



(n > 0) 



(n < 0) 



(70) 



with the domains of definition 



V 



s> 



X 



< V2m - In (1 + Ax) < - I 



(71) 



and V s< = [—1/A, 00) respectively. The corresponding phase trajectories are ruled by the 
expressions 





f 

1 






A 


exp | 


x(t) = < 








1 






A 


exp | 



A 



ay/2m 
A 



arcsm 



S-e 



£ 



cos + 00 



(H > 0) 



arcsinh 



£ 



-£ 



cos ( 2y/—£at + 



(H < 0) 



and 



p(t) 



\j2m,Q {£ — e) sin {2\/~£at + (, 
1 + Ax(t)) cos [V2^f In (1 + Ax(t))] 

^/2m (£ — e) sin {2\/—£at + O ) 



(n > 0) 



(H < 0) 



(72) 



(73) 



(1 + Ax(t)) cosh [v 7 ^! In (1 + Xx(t))] 

The potentials V s (x) and phase trajectories (x(t),p(t)) have been depicted in Figures [H 
and [2] for different values of the parameters. The description of the behavior of p(t) and 
m(x) is similar to that given in Section |3~T1 



3.3 Regular mass-functions 

The mass function 

m R {x) 



m 



1 + (Ax) 2 



(74) 
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V(x) V(x) v(.t) 




Figure 2: Upper line. Position dependent-mass hyperbolic Poschl- Teller potentials for 
the mass- functions m ss (left), m s (center), (right), and the parameters a = 1/2, 
= —1, mo = A = 1. Lower line. Some of the related (confined motion) classical 
trajectories in phase space for <p = 0, e = —2 and £ = — 1.5 (red curve), £ = — 1 (dashed 
curve), £ = —0.5 (blue curve). 



appears in the study of diverse oscillators including the nonlinear one [8,9J, and the related 
position-dependent mass versions [UEIIS] ■ This is a regular function defined on the whole 
real line. The Poschl- Teller potentials 



V R (x) 



cos 2 [\/2mof arcsinhAx] 



(n > o) 



(n < o) 



(75) 



are defined in 



V 



R> 



cosh 2 [y/2m f arcsinhAx] 

< v2m — arcsinhAa; < — > , 
A 2 J 



7T 

2 



(76) 
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and Vr k 



x(t) 



respectively. The phase trajectories read as 
A 



— sinh 
A 



sinh 



A 

Q\/2mo 



arcsm 



arcsinh 



£-6 

S 



cos ( 2y/Sat 



£-€ 



-£ 



cos ( 2y— Sat + 



(H > 0) 
(U > 0) 



(77) 



and 



Pit) 



'2m (S-e) sm(2VSat + 4> 
1 + (Xx(t)) 2 cos [f A/2m^arcsinh (Ax(t))] 

) sin (2 V /Z ^at + 0o) 



2m {£ 



1 + (Ax(t)) 2 cosh [f x /2m^arcsinh (Aar(*))] 



(ft > 0) 



(n < o) 



(78) 



The potentials Vr(x) and phase trajectories (x(t),p(t)) have been depicted in Figures [T] 
and [2] for different values of the parameters. The description of the behavior of p(t) and 
m(x) is similar to that given in Section [3TTT 

3.4 Exponential mass-functions 

The exponential mass-function m e (x), introduced in Eq. (|26|) . is a regular function defined 
on the whole real line. Here, we shall take k < in order to get finite masses in the positive 
regime of the domain. The domains of definition of the Poschl- Teller potentials 



VJx) 



cos 2 [y/2a 2 m Q ± (e KX / 4 - 1 



cosh 2 [V2a 2 m ^ (e 



kx/4 



{U > 0) 



(H < 0) 



(79) 



are respectively given by 

= [xq , x ], 



e> 



for any k < 0, m >0,l9a^0 



^e> = i x o j with 8y / 2c?mo < | «;| vr 



(80) 



and 2} 



e< 



with 



± 



ac 



IT 



7T AC 



8v2o^o 



Potential V e (x) in ( 179"]) is singular at for the domain T>\ y . This is also singular 

at and acquires a finite value as x — >■ +oo in the domain Pg>. Finally, this 
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potential cancels as x — > — oo, and goes to a finite value as x — > +00 in the domain D e< 
(see Figure [3]). The related phase trajectories are constructed according to 



x(t) = < 



' 4 






In 


K 

< 




4 






In 


K 





arcsin ( \ / - - cos(0 o + 2aV£ t ) + 1 



4 v / 2a 2 m 



4 v / 2a 2 m 



arcsinh 



ie-e 



cos 



6 + 2aV-£t + 1 



(ft > 0) 



(H < 0) 



(82) 



and 



p(t) 



e)m e Ka://4 sin( 



+ 2aV£t) 



cos [V2a 2 m ^ (e rea: / 4 - 1 



2{S - e)m e KX / 4 sin(0 o + 2a^£l) 
cosh [v / 2a 2 m 4 ; (e 



kx/4 



{n > 0) 



(^ < 0) 



The potentials V e (x) and phase trajectories (x(t),p(t)) have been depicted in Figure |3]for 
different values of the parameters. In this case we can distinguish three general kinds of 
potential: I. Potential V e (x) is such that scattering and positive bounded states (confined 
motion) are allowed (Figure El left). II. The potential is such that only confined motion is 
allowed (Figure [31 center). III. Negative and positive energies can be associated to scat- 
tering states, and negative bounded states are allowed (Figure [31 right). The description 
of the behavior of p(t) and m(x) is similar to that given in Section I3TT1 With the trajec- 
tories in the phase space (182]) and fl83|) we provide an explicit solution to the canonical 
problem involving "forces quadratic in the velocity" discussed in Refs. [TS] and 



4 Conclusions 

We have constructed the Lagrangian for a particle suffering a spatial variation of its mass. 
The corresponding Euler-Lagrange equations have been shown to recover the Newton's 
dynamical law associated to this system. As a consequence of the position-dependence of 
the mass, there is a force quadratic in the velocity which is non-inertial and represents 
the thrust of the system. The Lagrangian has been shown to be of the standard form 

2 

L = 2m(x) ~ "^( X )> m correspondence with the conditions studied in [TU]. In this form, the 

2 

construction of the related Hamiltonian also leads to the standard form W = 2 m{x) + V( X )- 
Thus, we have found that a simple description of this system consists of replacing the (con- 
stant) mass m by the appropriate function of the position m(x) in L and H. It has been 
shown that the canonical equations of motion must include a term including the thrust 
in order to recover the Newton's equation. This last result allows the identification of 
a canonical transformation mapping the variable mass problem to the one of a constant 
mass. Departing from the factorization of the Hamiltonian H = A + A~ + e, we have ob- 
tained two time-dependent integrals of motion Q which, in turn, allow the construction 
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Figure 3: Upper line. Position dependent-mass Poschl- Teller potentials for the expo- 
nential mass-function m e (x) in the domains of definition T> x e> (left), T> 2 e> (center), and 
T> e < (right). The parameters are a = 1, mo = 1/2, and e = 2, k = —2 (left), e = 2, 
k = —4 (center), e = —2, n = —4 (right). Lower line. Some of the related (confined 
motion) classical trajectories in phase space for <p = and £ = 2.5,4,5.5 (left and cen- 
ter), £ = —1.75, —1.5, —1 (right). In all cases red and blue curves correspond to reported 
lowest and highest energies. 



of the trajectories in the phase space (x(t),p(t)). Such invariants are associated to poten- 
tial functions which are necessarily of the Poschl- Teller form if one looks for the related 
spectrum generating algebras. The latter are closed by the Poisson structure defined in 
the factorizing of the Hamiltonian by demanding that the Poisson brackets of A and H. 
be expressed in terms of the polynomials of "H. Two general forms of the Poschl- Teller 
potentials have been found, one is of the trigonometric type and is connected with the 
su(l, 1) Lie algebra, the other is of hyperbolic type and is associated to the su{2) algebra. 
Finally, different solutions for the canonical equations of motion have been explicitly given 
in terms of particular forms of the mass-function m(x). We recover (as particular cases) 
Lagrangians and Hamiltonians already reported by other authors as associated to constant 
mass systems with forces quadratic in the velocity [T8]ll~9"]. Moreover, in contradistinction 
with such works, we provide explicit solutions for the corresponding equations of motion. 
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